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Abstract 

We report measurements of B meson decays to two kaons using 253 fb _1 of data collected with 
the Belle detector at the KEKB energy-asymmetric e + e~ collider. We find evidence for signals in 
B + — ► K°K + and B° — > K°K° with significances of 3.0a and 3.5a, respectively. (Charge-conjugate 
modes included) The corresponding branching fractions are measured to be B(B + — > K K + ) = 
(1.0 ± 0.4 ± 0.1) x 10~ 6 and B(B° i^°Z°) = (0.8 ± 0.3 ± 0.1) x 10" 6 . These decay modes are 
examples of hadronic b —> d transitions. No signal is observed in the decay B° — > K + K~ and we 
set an upper limit of 3.7 x 10~ 7 at 90% confidence level. 

PACS numbers: 11.30.Er, 12.15.Hh, 13.25.Hw, 14.40.Nd 
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Recent precise measurements of the branching fractions [l| and partial rate asymmetries 
[2] from the decays B — > Kn, txtx provide essential information to understand the B decay 
mechanism and to probe possible contributions from new physics. The rates for these decays 
constrain the hadronic b — > s and b — > u amplitudes. Here we report results on B° — > K°K 
and B + — > ~K°K + decays, which are examples of b — > d hadronic transitions. We also discuss 
a search for B° — > K + K~, which is sensitive to effects of final-state interactions (FSI) 0]. 
The results are based on a sample of 275 million BB pairs collected with the Belle detector 
at the KEKB e + e~ asymmetric-energy (3.5 on 8 GeV) collider |4| operating at the T(45) 
resonance. 

The Belle detector is a large-solid-angle magnetic spectrometer that consists of a silicon 
vertex detector (SVD), a 50-layer central drift chamber (CDC), an array of aerogel threshold 
Cherenkov counters (ACC), a barrel-like arrangement of time-of-flight scintillation counters 
(TOF), and an electromagnetic calorimeter (ECL) comprised of CsI(Tl) crystals located 
inside a superconducting solenoid coil that provides a 1.5 T magnetic field. An iron flux- 
return located outside the coil is instrumented to detect K\ mesons and to identify muons 
(KLM). The detector is described in detail elsewhere Two different inner detector 
configurations were used. For the first sample of 152 million BB pairs (Set I), a 2.0 cm 
radius beampipe and a 3-layer silicon vertex detector were used; for the latter 123 million 
BB pairs (Set II), a 1.5 cm radius beampipe, a 4-layer silicon detector and a small-cell inner 
drift chamber were used 6]. 

Charged kaons are required to have a distance of closest approach to the interaction point 
(IP) in the beam direction {£) of less than 4 cm and less than 0.1 cm in the transverse plane. 
Charged kaons and pions are identified using dE/dx information and Cherenkov light yields 
in the ACC. The dE/dx and ACC information are combined to form a K-tc likelihood ratio, 
1Z(K/ir) = Ck/{Ck + where Ck (At) is the likelihood that the track is a kaon (pion). 
Charged tracks with 1Z(K/tt) > 0.6 are regarded as kaons. Furthermore, charged tracks 
that are positively identified as electrons or muons are rejected. The electron identification 
uses information composed of E/p and dE/dx, shower shape, track matching % 2 , and ACC 
light yields, while information from the KLM, dE/dx and ACC are combined to identify 
muons. The kaon identification efficiency and misidentification rate are determined from a 
sample of kinematically identified D* + — > D°tc + , D° — » K~tt + decays, where the kaons from 
the D decay are selected in the same kinematic region as in B — > KK decays. The kaon 
efficiency is measured to be (84.24 ± 0.13)% for Set I and (82.84 ± 0.14)% for Set II, while 
the pion-fake-kaon rates are (5.40 ± 0.08)% and (6.86 ± 0.11)%, respectively. 

Candidate K° mesons are reconstructed through the Kg — > 7r + 7r _ decay. We pair 
oppositely-charged tracks assuming the pion hypothesis and require the invariant mass of 
the pair to be within 18 MeV/c 2 of the nominal Kg mass. Furthermore, the intersection 
point of the 7r + 7r _ pair must be displaced from the IP. 

Two variables are used to identify B candidates: the beam-constrained mass, Mt, c = 
J -EbLm — Pbi an d the energy difference, AE = E* B — E£ eaia , where E^ earn is the run depen- 
dent beam energy and E* B and p* B are the reconstructed energy and momentum of the B 
candidates in the center-of-mass (CM) frame, respectively. Events with Mb c > 5.20 GeV/c 2 
and \AE\ < 0.3 GeV are selected for analysis. 

The dominant background is from e + e~ — > qq (q = u, d, s, c) continuum events. Event 
topology and B flavor tagging information are used to distinguish between the spherically 
distributed BB events and the jet-like continuum backgrounds. We combine a set of modified 
Fox- Wolfram moments [^] into a Fisher discriminant. A signal/background likelihood is 
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formed, based on a GEANT-based |8[ Monte Carlo (MC) simulation, from the product of 
the probability density function (PDF) for the Fisher discriminant and that for the cosine of 
the angle between the B flight direction and the positron beam. The continuum suppression 
is achieved by applying a requirement on a likelihood ratio 71 = C S /(C S + £ gq -), where 
C s {qq) is the signal (qq) likelihood. Additional background discrimination is provided by B 
flavor tagging. For each event, the standard Belle flavor tagging algorithm provides 
a discrete variable indicating the probable flavor of the tagging B meson, and a quality 
r, a continuous variable ranging from zero for no flavor tagging information to unity for 
unambiguous flavor assignment. An event with a high value of r (typically containing a 
high-momentum lepton) is more likely to be a BB event, and a looser 1Z requirement can 
be applied. We divide the data into r > 0.5 and r < 0.5 regions. A selection requirement 
on 1Z for events in each r region of Set I and Set II is applied according to a figure-of-merit 
defined as iV s cxp / yj Ns xp + iVq| p , where iVJ xp denotes the expected signal yields based on MC 
simulation and the assumed branching fractions, 1.0 x 10~ 6 , and iV|| p denotes the expected 
qq yields from sideband data (Mb c < 5.26 GeV/c 2 ). 

Background contributions from T(4S) — > BB events are investigated using a large MC 
sample, which includes events from b — > c transitions and charmless decays. After all the 
selection requirements, no BB background is found for the B° — > K°K mode. Owing to 
K-tt misidentification, large B° — > K + tt~ and B + — > K°tt + feed-across backgrounds appear 
in the B° — ■> K + K~ and B + — > K°K + modes, respectively. A small charmless three-body 
contribution is found at low AE values for these two modes. 

The signal yields are extracted by performing unbinned two dimensional maximum like- 
lihood (ML) fits to the (Mb c , AE) distributions. The likelihood for each mode is defined 
as 

£ = exp (- N s,k,j)Y[(J2 N s,k,jP s ,k,j,i) 

s,k,j i s,k,j 

^s,k,j,i — Ps,k,j(M hci , AEi), 

where s indicates Set I or Set II, k distinguishes between events in the r < 0.5 and r > 0.5 
regions, i is the identifier of the z-th event, P(M^ C , AE) is the two-dimensional PDF of M^ c 
and AE, Nj is the number of events for the category j, which corresponds to either signal, qq 
continuum, a feed-across due to K-n misidentification, or background from other charmless 
three-body B decays. 

All the signal PDFs (P s ,fej=signai(Mb c , AE)) are parametrized by a product of a single 
Gaussian for M^ c and a double Gaussian for AE using MC simulations based on the Set 
I and Set II detector configurations. The same signal PDFs are used for events in the two 
different r regions. Since the Mb c signal distribution is dominated by the beam energy 
spread, we use the signal peak positions and resolutions obtained from B + — ► D°tt + data 
(TP — > Kg7i + 7T^ sub-decay is used for the K°K mode, while D° — > K + tt~ is used for 
the other two modes) with small mode dependent correlations obtained from MC. The 
MC-predicted AE resolutions are verified using the invariant mass distributions of high 
momentum D mesons. The decay mode D° — * K + it~ is used for B° — > K + K~ , D + — > Kg7r + 
for B + — ► K°tt + and D° — > KgTT + 7T^ for B° — > K°K° . The parameters that describe the 
shapes of the PDFs are fixed in all of the fits. 

The continuum background in AE is described by a linear function while the M^ c distri- 
bution is parameterized by an ARGUS function f(x) = xy/1 — x 2 exp [— £(1 — x 2 )], where 
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x is Mb c divided by half of the total center of mass energy Therefore, the continuum 
PDF is the product of this ARGUS function and the linear function, where the overall nor- 
malization, £ and the slope of the linear function are free parameters in the fit. These free 
parameters are r-dependent and allowed to be different in Set I and Set II. The background 
PDFs for charmless three-body B decays for the K + K~ and K K + modes are each modeled 
by a smoothed two-dimensional histogram, obtained from a large MC sample. The feed- 
across backgrounds for these two modes from the K + tt~ and K°7r + events have M^ c — AE 
shapes similar to the signals with the AE peak positions shifted by ~ 45 MeV. The methods 
to model the K + K~ and K K + signal PDFs are also applied to describe the feed-across 
background. 

When likelihood fits are performed, the yield for each background component {N S} kj 
where j = qq, feed-across, charmless) is allowed to float independently for each s (Set I or 
Set II), and k bin (low or high r region). For the signal component, the same branching 
fraction is required by constraining the number of signal events in each (s, k) bin using the 
measured efficiency in the corresponding (s, k) bin. Table U summarizes the fit results for 
each mode. We observe 13.3 ± 5.6 ± 0.6 K°K + and 15.6 ± 5.8lJ;g K°TC° signal events with 
significances of 3.0<r and 3.5a, respectively. The second errors in the yields are the systematic 
errors from fitting, estimated from the deviations after varying each parameter of the signal 
PDFs by one standard deviation, and from modeling the three-body background, studied by 
excluding the low AE region (< —0.15 GeV) and repeating the fit. At each step, the yield 
deviation is added in quadrature to provide the fitting systematic errors and the statistical 
significance is computed by taking the square root of the difference between the value of 
— 21n£ for the best fit value and zero signal yield. The smallest value is chosen to be the 
significance including the systematic uncertainty. 

Figure Q shows the M^ c and AE projections of the fits after requiring events to have 
\AE\ < 0.06 GeV and 5.271 GeV/c 2 < M bc < 5.289 GeV/c 2 , respectively. The feed- 
across yields are 47.1 ± 8.7 in the K + K~ mode and 16.4 ± 6.1 in the K°K + mode. The 
amounts of the feed-across background are consistent with the expectations of 49.1 K + 7r~ 
and 18.8 K°n + events, based on MC simulation and measured branching fractions [l2^ . 
The MC modeling of the requirement on the likelihood ratio, 1Z is investigated using the 
B + -> D°7i + (D° -> K° s n + 7r- for K°K° and D° -> K+ir~ for the others) samples. The 
obtained systematic errors are ±2.9% for B° — > K°K° and ±6.8% for the other two modes. 
The systematic error on the charged track reconstruction efficiency is estimated to be around 
1% per track using partially reconstructed D* events. The resulting Kg reconstruction is 
verified by comparing the ratio of D + — > K s n + and D + — > K~7t + 7t + yields with the MC 
expectation. The resulting K s detection systematic error is ±4.5%. The final systematic 
errors are then obtained by quadratically summing the errors due to the reconstruction 
efficiency and the fitting systematics. 

With 275 million BB pairs, we find evidence of B+ -> K K + and B° -> K°K with 
branching fractions B(B + -> K°K + ) = (1.0 ± 0.4 ± 0.1) x 10~ 6 and B(B° -> K°K°) = 
(0.8±0.3±0.1) x 10~ 6 . These are examples of hadronic b — » d transitions. Our measurements 
are consistent with preliminary results reported by the BaBar collaboration and agree with 
some theoretical predictions [13, Hi llR . It has been suggested that the branching 
fraction and CP asymmetry of the mode B° — > K°~K° , which originates from the flavor- 
changing neutral current process b — > dss, may be sensitive to physics beyond the Standard 
Model [l6j. Measurements with larger statistics are needed for this purpose. No signal is 
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FIG. 1: M hc (left) and AE (right) distributions for B° K+K~ (top) and 5+ -> K K + 
(middle) and i?° — > K°K candidates. The histograms show the data, while the curves represent 
the various components from the fit: signal (dashed), continuum (dotted), three-body B decays 
(hatched), background from mis-identification (dash-dotted), and sum of all components (solid). 
In the K + K~ mode, there is a large contribution from misidentified K + tt~ but no significant signal 
excess. In the K°K + mode, the signal and misidentified K°ir + contributions are comparable in 
size. In the K°K° mode, there is a signal excess but no misidentification background. 



observed in B° — ► K + K~ and we set the upper limit of 3.7 x 10~ 7 at the 90% confidence 
level, using the Feldman-Cousins approach taking into account both the statistical and 
systematic errors jl9^ . 

We thank the KEKB group for the excellent operation of the accelerator, the KEK 
cryogenics group for the efficient operation of the solenoid, and the KEK computer group 
and the Nil for valuable computing and Super-SINET network support. We acknowledge 
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TABLE I: Fitted signal yields, reconstruction efficiencies, product of efficiencies and sub-decay 
branching fractions (B s ), branching fractions and significances for individual modes. 



Mode 


Yield 


Eff.(%) ES.xB s ( 


%) £(10- 6 ) 


Sig. 


K+K- 


2 5+ 5 - 1 


15.5 15.5 


< 0.37 


0.5 


K°K+ 


13.3 ±5.6 


14.5 5.0 


1.0 ±0.4 ±0.1 


3.0 


K°K° 


15.6 ±5.8 


28.7 6.8 


0.8 ±0.3 ±0.1 


3.5 
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